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Two paradigm shifts in DNA sequencing technologies — fi-om bulk to single molecules and from optical to 
electrical detection — are expected to reahze label-free, low-cost DNA sequencing that does not require PGR 
amplification. It will lead to development of high- throughput third-generation sequencing technologies for 
personalized medicine. Although nanopore devices have been proposed as third-generation 
DNA-sequencing devices, a significant milestone in these technologies has been attained by demonstrating a 
novel technique for resequencing DNA using electrical signals. Here we report single-molecule electrical 
resequencing of DNA and RNA using a hybrid method of identifying single-base molecules via tunneling 
currents and random sequencing. Our method reads sequences of nine types of DNA oligomers. The 
complete sequence of 5'-UGAGGUA-3' from the let-7 microRNA family was also identified by creating a 
composite of overlapping fragment sequences, which was randomly determined using tunnehng current 
conducted by single-base molecules as they passed between a pair of nanoelectrodes. 

The Human Genome Project has provided a complete genetic road map, which has expanded our knowledge 
about not only biology of humans but also of other organisms. We believe that the emerging high- 
throughput DNA sequencing technologies would allow us to acquire such information at a much greater 
rate than those possible using the current technologies. For example, we will be able to follow up on the evolution 
of viral and bacterial drug resistance in real time, uncover the huge diversity of novel genes that are currently 
inaccessible, and provide physicians with genomics-based preventive diagnostics and therapeutics that specif- 
ically and comprehensively target individual patients' \ 

First and second generation DNA sequencing technologies identify nucleotides through light emission. They 
require PGR for the amplification of sequencing templates in order to produce sufficient material to generate 
detectable signals. Third generation DNA sequencing technologies detect single nucleotides directly by changes 
in electric current, thereby negating the need for either fluorescent probes or PGR amplifications' ''. 

Nanopores are a key component of third generation DNA sequencing technologies'" These devices provide two 
possible approaches to DNA sequencing. The first includes use of bio-nanopores^ " or solid-state nanopores'". As 
each nucleotide in a DNA strand translocates through a nanopore, it specifically modulates the ionic current being 
passed through the nanodevice; this allows the identification of each nucleotide, and eventually of the DNA 
sequence"""'. The second method includes use of gating nanopores'""^'', in which DNA can be identified via the 
tunnelling current they conduct as they pass between a pair of nanoelectrodes^*"^'' (Fig. la). Measurable changes in 
the current occur in response to small differences in the electronic structure of the DNA molecule, making it 
possible to identify the type and arrangement of individual base molecules. 

Single-molecule electrical sequencing via tunnelling current provides two means of sequencing long, single- 
stranded DNA molecules. One method enables one-by-one identification of the base molecules along a long, 
single-stranded DNA molecule"'^. The alternative method enables the determination of the sequence of a DNA 
molecule as a composite of overlapping randomly-fragmented sequences; this method is aptly termed as random 
DNA sequencings''^". 

Because single-molecule electrical-random sequencing has not yet been reported, the determination of random 
sequences of short DNA and RNA composed of known sequences (resequencing DNA and RNA) is a significant 
milestone in sequencing. Here we report about single-molecule electrical resequencing of DNA and RNA by a 
hybrid method of identifying single base molecules via a tunnelling current and random sequencing. Nine DNA 
oligomers composed of 3 base molecules were sequenced by our single-molecule electrical sequencing method. 
5'-UGAGGUA-3' miRNA was completely resequenced by our hybrid method of single-molecule electrical 
sequencing and random resequencing. 
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Figure 1 | Schematic diagram of single-molecule electrical sequencing of 
DNA and RNA. (a) Single base molecules on a DNA molecule in a solution 
are identified via a tunnelling current that passes between nanoelectrodes. 
Nanogaps of 0.8 nm were fabricated using a nano-MCBJ. A constant DC 
voltage of 0.4 V was applied between the nanoelectrodes, and electrical 
signals were measured at a sampling rate of 10 KHz. (b) Molecular 
structures of deoxynucleoside monophosphates of DNA and (c) 
nucleoside monophosphates of RNA. 



Results 

Single-molecule identification of DNA and RNA bases. In this 
study, we investigated the conductance values of the four bases of 
DNA and RNA. Four deoxyribonucleoside monophosphates (dAMP, 
dCMP, dGMP, dTMP) and four ribonucleoside monophosphates 
(rAMP, rCMP, rGMP, rUMP) could be identified by the tunnelling 
current they conducted as they passed between a pair of nanoelec- 
trodes separated by a gap of 0.8 nm (Fig. 1). This gap was set using 
a nanofabricated, mechanically controllable break junction (nano- 
MCBJ)""' (Supplementary Figs. S1-S6). We also determined the time 
dependence of the electrical current between the nanoelectrodes when 
each of the single nucleotides stochastically passed through the nano- 
electrode gap by Brownian motion. Spike-hke electrical signals were 
observed at high frequencies for approximately 100 ms and peak ele- 
ctrical currents ranged from 10 pA to several hundred pA, indicating 
the presence of a base molecule between the nanoelectrodes (Figs. 2a 
and 2b and Supplementary Fig. S7). The electrical signals were 
characterised by Ip and f^, which were defined as the maximum 
current and the duration of the current, respectively. Typical values 
of 7p and for dGMP were 100 pA and 1 ms, respectively (Sup- 
plementary Fig. S8). 

Single-molecule conductance (/p/V) histograms were then con- 
structed from the Ip data using approximately 1000 spikes for each 
molecule and were analysed using Gaussian fit curves. The G values 
(the peak value in each conductance histogram) for each single nuc- 
leotide was found to be in the following order: dGMP (87 pS) > 
dAMP (67 pS) > dCMP (60 pS) > dTMP (39 pS) for DNA and 
rGMP (123 pS) > rAMP (92 pS) > rCMP (64 pS) > rUMP (50 pS) 
for RNA. The relative G normalized by dGMP and rGMP was also 
compared to explore the order of conductance and the following 
results were found: dGMP = rGMP = 1, dAMP (0.77) ~ rAMP 
(0.75) > dCMP (0.69) > rCMP (0.52) > dTMP (0.45) > rUMP 
(0.41) (Table 1). Calculations based on density functional theory 
indicated that the order based on the highest occupied molecular 
orbital (HOMO) energy was as follows: guanine ( — 5.7 eV) > aden- 
ine ( — 5.9 eV) > cytosine ( — 6.1 eV) > thymine ( — 6.6 eV) > uracO 
( — 6.9 eV) (Supplementary Table SI). Note that this order corre- 
sponds to the order of the relative G values, suggesting that our 
single-molecule electrical detection method can identify molecular 
species based on characteristic energy levels, particularly the HOMO 
energy level. 

Single-molecule electrical resequencing of DNA oligomers. In 

addition, dual level conductance-time profiles were obtained for 
TGT, GTG, ATA, GAG and GAG DNA oligomers (Figs. 3b-3f and 
Supplementary Fig. S9). Single-molecule conductance histograms 
constructed from the G-f profiles showed two peaks corresponding 
to dual level conductance. When the relative G for the higher G 
corresponding to single-molecule conductance of dGMP was set to 
1, the relative G for the lower G in TGT, GTG, and GAG conductance 
histograms was calculated as 0.29 ± 0.12, 0.35 ± 0.12, and 0.68 ± 
0.08, respectively. (Table 2). Since these values were within the ranges 
of the relative G for dTMP (0.45 ± 0.12) and dAMP (0.77 ± 0.20), 
the obtained values, 0.29, 0.35 and 0.68, were assigned to dTMP, 
dTMP and dAMP, respectively. Similarly, when the higher G was 
set relative to G = 0.77 (obtained from the conductance histogram 
for dAMP), the lower relative G in ATA and CAC conductance 
histograms was calculated to be 0.41 ± 0.07 and 0.52 ± 0.11, 
respectively. Since these values were within the ranges of the 
relative G for dTMP (0.45 ± 0.12) and dCMP (0.69 ± 0.25), the 
obtained values, 0.41 and 0.52, were determined to be those of dTMP 
and dCMP, respectively. Therefore, the type of base molecules 
forming DNA oligomers could be identified using the relative G of 
the single base molecules. 

To resequence the base molecules of a single DNA or RNA mole- 
cule, we had to extract information about the arrangement of the 
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Figure 2 | Single-molecule identification of deoxynucleoside monophosphates of DNA and nucleoside monophosphates of RNA. (a) Typical current- 
time profiles obtained from 0.10 |iM solutions of dGMP. (b) Expanded figure of a typical electrical signal, which is characterized by and 7p, where /p is 
defined as the maximum current and as the duration of the current, (c) Conductance histograms of four deoxynucleoside monophosphates of DNA, 
constructed from 1000 7p for each molecule. The blue, red, purple and green lines indicate the single-molecule conductance of dGMP, dAMP, dCMP and 
dTMP, respectively. The colour codes of lines corresponding to the Jp histograms represent Gaussian fitted lines, (d) Conductance histograms of four 
nucleoside monophosphates of RNA. The blue, red, purple and green lines indicate the single-molecule conductance values of rGMP, rAMP, rCMP and 
rUMP, respectively. The colour codes of lines corresponding to the 7p histograms represent Gaussian fitted lines. The single-molecule conductance values 
are summarized in Table 1. The 7p histograms were constructed from 1000 Ip data points for each nucleoside monophosphate, (e) Typical t<j histograms of 
dGMP and GGG constructed from 1000 data. The line curve shows a Gaussian fitted line with a peak tj of approximately 1 ms. 
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Table 1 | Single-molecule conductance (G) and relative single- 
molecule conductance of deoxynucleoside monophosphates 
(dGMP, dAMP, dCMP and dTMP) of DNA and nucleoside mono- 
phosphates (rCMP, rAMP, rCMP rUMP) of RNA. Single-molecule 
conductance were obtained from the conductance histograms in 
Figs. 2c and 2d. Single-molecule conductance was normalized to 
that of dGMP and rGMP. FWHM, full width at half maximum 
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bases from the /- 1 profiles. Since the electric signals were composed of 
/p and tjj and the type of base molecule was determined from /p, 
the base arrangement could be identified using t^. Here we investi- 
gated the I-t profiles of dGMP and the GGG DNA oligomer 
(Supplementary Fig. SIO. Supplementary Figure Sll shows I-t pro- 
files of AAA, CCC, and TTT). One-level spike-liice signals were 
observed in the I-t profiles of GGG as well as dGMP (Fig. 3a). The 

histogram of dGMP showed a peak at approximately 0.8 ms, while 
that of GGG gave a peak at approximately 0.8 ms and extended to 
1000 ms (Fig. 2e and Supplementary Fig. SB). Therefore, it was clear 
that extraction of electrical signals with a t^ of more than 2 ms 
(—0.8 ms X 3) was required. Using this criterion, electrical signals 
of single DNA molecules consisting of n base molecules should show 
n plateaus. For example, the electrical signals of G and T, GT and TG, 
TGT and GTG would be expected to show one, two and three pla- 
teaus, respectively. Thus, the criteria for identifying electrical signals 
from single DNA molecules consisting of three base molecules are 
that td should be >2 ms and that the number of plateaus should be 3. 

Electrical signals (Figs. 3b-3f) were then automatically extracted 
using the criteria described above and three obvious plateaus were 
observed (Supplementary Figs. S12 and S13, Tables S2 and S3) . In 
one case, two lower plateaus of T were observed in the first and third 
positions, with G forming the middle higher plateau (TGT) (Fig. 3g); 
in the other case, G was observed to form the first and third higher 
plateaus, with T in the lower centre position (GTG) (Fig. 3h). The 
results clearly demonstrated that DNA oligomers could be rese- 
quenced. In the experimental conditions using the nano-MCBJ not 
only GTG and TGT but also G, T, TG, GT, GTGTT and TGTGT 
could be identified because single DNA molecules were stochastically 
trapped between the nanoelectrodes because of Brownian motion. 
For example, when Brownian motion turns over a single TGT oli- 
gomer at the T position between nanoelectrodes, we can obtain 
unexpected G-t profiles such as TGTGT composed of a forward 
sequence (TGT) and a backward sequence (GT). Therefore, single- 
molecule electrical resequencing using nano-MCBJ can be used to 
randomly identify sequences of single base DNA molecules. 

Single-molecule electrical random resequencing of miRNA. In 

random sequencing methods^^'^", a DNA molecule is broken up 
randomly into numerous small fragments. Multiple overlapping 
reads for the target DNA are obtained by performing several 
rounds of fragmentation and sequencing. Computer programs 
then use the overlapping ends of different randomly fragmented 
sequences to assemble them into continuous sequences (sequence 
contigs) and further to assemble the contigs into a continuous 
sequence. To assist in the effort to evaluate the potential for single- 
molecule electrical resequencing, we selected a common sequence 
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Figure 3 | Resequencing of DNA oligomers. Relative conductance 
(G)-time profiles of (a) GGG, (b) TGT, (c) GTG, (d) ATA, (e) CAC and 
(f) GAG. The relative G values are summarized in Table 2. The left panels in 
each figure show the related G histograms. Upper and lower panels show 
expanded 7- 1 and smoothed G-tprofiles of (g) TGT and (h) GTG. The blue, 
red, purple and green lines and bands in all figures represent the relative G 
and its dGMP, dAMP, dCMP and dTMP bands, respectively. 

(5'-UGAGGUA-3') from the let-7 miRNA family^", which is a well- 
known cancer marker and consists mainly of 22 bases. 

First, when we measured the I-t profiles of miRNA, the current 
histogram constructed from the profile showed three peaks at / = 70 
pA, 50 pA, and 33 pA (Fig. 4a). The peaks were assigned to rGMP, 
rAMP and rUMP, since the relative G values of 1, 0.71 and 0.47 agree 
with those of these single base molecules for RNA (Table 3). This 
result shows that three base molecules exist in the miRNA. Fragment 
sequences were determined using the same procedure as was used to 
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Table 2 | Relative single-molecule conductance of deoxyribonucleotides obtained from current-ti 
FWHM. Reference shows the relative single-molecule conductance extracted from Table 1 
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dGMP 1 .00 ± 0.25 
dAMP 0.77 ± 0.20 
dCMP 0.69 ± 0.25 
dTMP 0.45 ±0.12 


1.00 ±0.29 1.00 ±0.18 
0.35 ± 0.12 0.29 ±0.12 


0.77 ±0.18 
0.41 ± 0.07 


0.77 ±0.14 
0.52 ±0.12 


1.00 ±0.19 
0.68 ± 0.12 



resequence DNA oligomers. In particular, electric signals fluctuating 
with a period of less than 0.8 ms were ignored because the character- 
istic dwell time for single bases was 0.8 ms, as shown in Fig. 2e. 
Figures 4b-4d show typical analysed electrical signals. As is clearly 
shown in the figures. A, G, U, AU, UGAGG and UGAGGUA were 
identified as fragment sequences (Supplementary Fig. SI 4). When 
133 electrical signals were analysed in the same way, we obtained 
randomly fragmented sequences, including 19XA, 15XG, 44XU, 
5XUA, lOXGA, 5XUG and the 35 fragments shown in Fig. 4e. As 



mentioned when resequencing DNA oligomers, characteristic frag- 
ment sequences (e.g. GAGAGGUA, UGAGGAGA and UGAGGU- 
AUA) were obtained because of stochastic motions originating from 
Brownian motion. In addition, both the forward and reverse reads 
were observed. For example, out of the 15 asymmetric reads of more 
than three bases (e.g. UGA and GAGGUA), the forward (UGA and 
GAGGUA) and reverse (AGU and GAGGUA) had yields of 60% 
and 40%, respectively. Furthermore, we acquired misread frag- 
ments of AGAUA and GAGGUG corresponding to AGGUA and 
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Figure 4 | Random resequencing of 5'-UGAGGUA-3' microRNA. (a)Typical J-f profiles from a 1.0 solution of UGAGGUA with 100 mM 
phosphate buffer, (b)-(d) Typical relative G histograms, raw and smoothed G- tprofiles of randomly fragmented sequences. Purple, red and blue lines and 
bands show relative G values (Table 3) and rUMP, rAMP and rGMP bands, respectively. (U, G, A, UA), (AU, UGAGG) and (U, UA, UGAGGUA) were 
resequenced according to the G-f profiles in (b), (c) and (d), respectively, (e) The 35 randomly fragmented sequences used for resequencing the 
microRNA. Sequence contigs 1-4 were assembled by fragment sequences within the red, green, blue and purple frames. The full sequence of the 
microRNA was determined by overlapping the 4 sequence contigs. 
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Table 3 | Relative single-molecule conductance of ribonucleotides 
in random fragments of 5'-UGAGGUA-3', shown in Figs. 4b-4d. 
Reference shows the relative single-molecule conductance 
extracted from Table 1 

Reference Fig. 4b Fig. 4c Fig. 4d 

rGMP 1.00 ±0.44 1.00 ±0.16 1.00 ±0.13 1.00 ±0.13 
rAMP 0.75 ±0.27 0.65 ± 0.15 0.55 ± 0.16 0.73 ±0.16 
rUMP 0.41 ±0.10 0.33 ±0.07 0.29 ±0.11 0.36 ± 0.09 



GAGGUA, respectively, because the relative G band of rGMP over- 
laps that of rAMP. 

Next, we used the overlapping ends of the different randomly 
fragmented sequences to assemble these into continuous sequences 
(sequence contigs) and then to assemble the contigs into a continu- 
ous sequence. One and two read signals (short signals) were too short 
to determine their assigned positions on the sequence of miRNA. On 
the other hand, the longer signals (over three reads) can be assigned 
on the sequence of the miRNA, as opposed to the known sequence of 
the miRNA. Therefore, only the longer signals were used for signal 
assembly and formation of contigs. The 35 fragment sequences were 
used to assemble the sequence contigs (Fig. 4e). Sequence contig 1 
was assembled as UGA using 13 fragment sequences within the red 
frame in Fig. 4e. Similarly, sequence contigs 2-4 were assembled as 
GAGG, AGGUA and AGGU, respectively, using 17, 10 and 13 frag- 
ment sequences within the green, blue and purple frames in the same 
figure. Since there was no sequence gaps between the sequence con- 
tigs due to the short miRNA used, the four sequence contigs were 
assembled as UG AGGUA. The read error was 10% because by using 
35 fragments, the number of wrong read bases (gray bases shown in 
Fig. 4e) was 17 out of a total of 168 read bases. 

Discussion 

Single base molecule can be identified using G in DNA and RNA 
monomers, though conductance histograms exhibit significant over- 
lap due to the large full width at half maximum (FWHM) of G 
(Figs. 2c and 2d). The large FWHM largely originates from the fact 
that tunneling current is highly sensitive to configurations of single 
base molecules with respect to the nanoelectrodes. Therefore, in an 
effort to decrease the FWHM and to improve the sequencing 
reliability, it is important to decrease the degree of freedom of 
configurations. 

The electrical measurement of DNA monomers, DNA oligomers 
and a miRNA showed interesting results to provide one method to 
reduce the FWHM. We found large FWHM of DNA and RNA 
monomers, where the FWHM of dGMP, dAMP and dCMP was 
more than 0.40 while that of rGMP and rAMP were 0.88 and 0.54, 
respectively (Table 1). However, the FWHM of dGMP in TGT and 
GAG oligomers decreased though that in GTG showed a slight 
increase (Table 2). The FWHM of dAMP, dCMP, dTMP in 4 
DNA oligomers reduced, and that of dCMP in CAC oligomer 
reduced to 48% of that of single dCMP molecule. Furthermore, frag- 
mented sequences of miRNA had smaller FWHM (Table 3). In par- 
ticular, the largest FWHM (0.88) of rGMP in the miRNA among 8 
base molecules reduced to 36%. The FWHM of rAMP in the miRNA 
also reduced to 59% of that in the single base molecule. As a result, 
conductance histograms of rGMP, rAMP and rUMP in fragmented 
sequences exhibited no overlap within the FWHMs (Figs. 4b-4d). 

Reduction of FWHM of base molecules in DNA oligomers and a 
miRNA comes from decreasing the degree of freedom of configura- 
tions, where chemical bonds via sugars and phosphates restrict con- 
formations of single base molecules. In this study, single molecules 
could move around nanoelectrodes fabricated by nano-MCBJ. 
Therefore, using gating nanopores composed of nanoelectrodes 
and nanopores, we would be able to obtain smaller FWHM because 



configurations of single DNA and RNA molecules are strongly 
restricted by nanopores. 

Methods 

Fabrication of nanogaps. A microfabricated MCBJ was employed to form 
nucleotide-sized electrode gaps (Supplementary Figs. SI and S2). The fabrication 
procedures for MCBJs are described elsewhere.'^^ First, we coated a phosphor bronze 
substrate with a thin polyimide layer for electrical insulation. A gold nanoj unction 
was then fabricated using electron beam lithography and a subsequent lift-off process. 
The junction was then exposed to isotropic reactive ion etching using CF4/O2 gas to 
remove the underlying polyimide and obtain a free-standing gold nanobridge. The 
MCBJ sample was mounted in a three-point bend configuration. The substrate was 
then bent and the j unction broken mechanically to form a pair of gold nanoelectrodes. 
The junction conductance was monitored using a Keithley 6487 picoammeter 
(Keithley) under a DC bias voltage of 0.1 V. A lO-kH resistor was connected in series 
during this process to protect the junction from overcurrent-induced 
electromigration failure. Special care was taken while adjusting the electrode gap size 
by a self-breaking technique. 

Materials. We purchased dAMP (2'-deoxyadenosine-5'-monophosphate: Sigma- 
Aldrich), dCMP (2'-deoxycytidine-5'-monophosphate sodium salt: Sigma- Aldrich), 
dGMP (2'-deoxyguanosine-5'-monophosphate sodium salt hydrate: Sigma- 
Aldrich), dTMP (Thymidylic acid disodium salt: Tokyo Chemical Industry Co. 
(TCI)), rAMP (2 '-adenosine-5' -monophosphate disodium salt: Oriental yeast), 
rCMP (cytidine 5 '-monophosphate disodium salt: TCI), rGMP (Guanosine 
5 '-monophosphate sodium salt hydrate: TCI) and rUMP (uridine 
5 '-monophosphate disodium salt hydrate: TCI). These were used to make 10 [iM 
aqueous solutions without further purification. DNA oligomers were synthesized by 
FASMAC Co., Ltd. and Hokkaido System Science Co., Ltd. RNA oligomers were 
synthesized by Operon Biotechnologies Co., Ltd. No buffer was added to the sample 
solutions of single nucleotides. A total of 100 mM phosphate buffer was used for all 
the sample DNA and RNA oligomer aqueous solutions. 

Test procedure. We formed a 0.8-nm electrode gap in a dilute Milli-Q solution of a 
target nucleotide molecule at a concentration of 0.10 [.iM. The current across the 
electrodes was recorded at 10 kHz using a custom-built logarithmic current amplifier 
and a PXI-4071 digital multimeter (National Instruments) under a DC bias voltage of 
0.4 V. After every 1 h of 7-t measurement, we replaced the MCBJ sample with a new 
one. With the single-nucleotide measurements, in order to suppress a possible 
contamination effect on the measurements, about 20 samples were used for each 
nucleotide to obtain approximately 1000 spike-like signals for constructing the /p 
histograms. For the oligonucleotide measurements, about 10 samples were used for 
each sample oligonucleotide to obtain approximately 200 signals. 

Resequencing procedure. Resequencing was performed using the following 
procedure: (1) measurement of 1-t profiles; (2) extraction of electrical signals with 
id > 2 ms; ( 3 ) plotting of 7 histogram using all data of 7- f profiles; (4) determination of 
current baseline on the basis of the lowest conductance peak; (5) smoothing of all data 
in an attempt to obtain plateaus; (6) determination of base-type molecules using the 
relative single-molecule conductance of base molecules, as shown in Table 1; (7) base 
assignment based on relative single-molecule conductance and plateaus. The detailed 
sequencing procedures for GGG and GTG are shown in Supplementary Figs. S 12 and 
S13, respectively. 
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